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ABSTRACT 


This  report  describes  a  computer  program  for  determining  the 
backscattered  fields  of  a  conducting  body  of  revolution  composed  of 
sections  of  cones  and  cylinders.  The  target  may  be  closed  at  one  or 
both  ends  with  circular  disks.  The  target  may  have  as  many  as  20 
sections,  and  the  program  can  readily  be  modified  to  handle  a  larger 
number.  The  backscattered  field  for  Eq  or  polarization  is  computed 
using  wedge  diffraction  theory  and  geometrical  optics.  The  computed 
results  are  in  good  agreement  with  experimental  measurements  for  cones, 
cylinders,  double  cones,  and  conically  capped  cylinders. 
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A  COMPUTER  PROGRAM  FOR  BACKSCATTER  BY  TARGETS 
COMPOSED  OF  CONES,  CYLINDERS,  AND  DISKS 


INTRODUCTION 

A  computer  program,  which  computes  the  backscattered  field  of  a 
body  of  revolution  composed  of  cones,  frustums,  cylinders,  and  disks 
is  given  in  Appendix  I.  This  program  is  coded  in  Fortran  IV  and  the 
cards  are  numbered  consecutively  in  columns  76-80.  This  program  has 
been  tested  for  cones,  cylinders,  frustums,  conically  capped  cylinders, 
and  double  cones,  for  Ey  polarization  of  the  incident  and  scattered  fields. 
The  results  of  these  tests  are  presented  in  Reference  1.  The  theoretical 
basis  for  this  program  is  described  in  References  1  and  2. 

The  function  of  this  program  is  to  compute  the  backscattered 
field  as  a  function  of  aspect  angle  and/or  frequency  for  a  given  target. 
This  is  accomplished  by  reading  and  storing  the  target  description, 
identifying  the  regions  of  aspect  angle  corresponding  to  axial  and 
specular  directions,  and  applying  the  appropriate  geometrical  theory  of 
diffraction  solutions  to  calculate  the  scattered  fields.’  In  the  form 
presented  here  the  angular  pattern  for  a  given  target  and  frequency  is 
computed.  Modification  of  the  program  to  compute  a  frequency  curve 
for  a  given  target  and  aspect  angle  is  straightforward. 


TARGET  DESCRIPTION 

The  target  shape  is  described  in  cylindrical  coordinates  by  the 
second  degree  equation 


( 1 )  F(p,  z)  =  Aip2  +  Azt}  +  Asp  z  +  A4Z  +  A5P+A6-O  • 


As  this  program  is  written  to  handle  only  targets  composed  of  cones, 
frustums,  cylinders  and  disks,  the  description  of  the  profile  of  the 
target  is  that  of  a  straight  line  segment,  and  may  be  expressed  as 


(2)  F(p,  z)  =  A4Z  +  A5  p  +  A6  -  u 

where  Ax  =  A2  =  A3  =  0  .  In  addition  to  Eq.  (2)  the  boundaries  of  each 
section  of  the  target  must  be  specified.  That  is,  within  each  section  of 
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the  target  the  profile  is  specified  by  a  set  of  constants  A4,  A5  ,  .  The 

next  section  is  described  by  a  new  set  of  constants  and  so  forth.  As  an 
example  consider  the  target  shown  in  Fig.  2.  This  target  is  composed 
of  a  cone,  a  cylinder,  and  a  disk,  and  is  a  three  section  target.  Thus 
a  set  of  constants  must  be  specified  for  each  section  and  the  angular 
boundaries  of  the  sections  specified.  Figure  1  gives  the  relations  by 
which  the  constants  may  be  determined  for  each  surface.  The  computer 
format  for  this  input  data  will  be  discussed  below. 

If  Eq.  (1)  is  used  to  describe  the  target,  i.e.  any  of  the  constants 
A\>  Az  A3  non- zero,  this  program  will  compute  the  scattered  field  due 
to  the  wedge  type  discontinuities  on  the  target  but  will  not  calculate  the 
correct  total  scattered  field.  This  is  because  there  is  no  provision  in 
this  program  to  calculate  the  geometrical  optics  field  of  a  doubly  curved 
surface.  A  calculation  of  this  geometrical  optics  field  is  included  in 
the  creeping  wave  computer  program.3 

WEDGE  DIFFRACTION  COMPUTER  PROGRAM 

Referring  to  the  computer  program  listing  shown  in  Appendix  I, 
the  function  of  the  significant  sections  of  the  program  will  be  discussed. 
The  card  numbers  associated  with  each  section  will  be  specified.  This 
discussion,  together  with  the  comment  cards  included  in  the  program 
listing,  is  intended  to  give  sufficient  information  about  the  program  to 
enable  a  qualified  programmer  to  both  use  and  modify  the  program. 
Statements  which  are  in  common  use  in  Fortran  IV  such  as  DIMENSION, 
COMPLEX,  and  FORMAT  statements  will  not  be  discussed  as  it  is 
assumed  that  the  reader  has  a  knowledge  of  Fortran  IV. 

The  COMMON  declaration  (0006)  is  used  to  store  the  constants 
required  in  Eq.  (2)  in  the  common  block  labelled  /  DATA/.  This  com¬ 
mon  block  is  used  in  conjunction  with  the  unlabelled  common  block  to 
transfer  a  particular  set  of  constants  Al(I)  to  A6(I)  into  the  unlabelled 
common  regions  shared  by  the  subroutines.  This  provision  reduces  the 
number  of  calling  variables  required  by  each  subroutine. 

The  READ  (0030-0040)  statements  in  this  block  of  statements  read 
the  required  data  and  provision  is  also  made  to  write  out  this  data  for 
the  purpose  of  identification. 

The  statements  0041-0045  initialize  constants  which  are  required 
in  the  calculations. 
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£  *  TAN  1  (  ±) 

b 

A4  -  ' 

a5  *  0 


£  =  180°  -  TAN  1  (-£-) 

D 
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A6  =  -b 


Fig.  1(a).  Straight  line  segments  and  their  corresponding 
analytic  description. 
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Fig.  1(b).  Straight  line  segments  and  their  corresponding 
analytic  description. 
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a  =  TAN  1  (y)  p  =  180°  -  TAN  1  ( -g-  ) 
A4  =  0  1 

A5  =  1 

a6  *  -a 

Fig.  1(c).  Straight  line  segments  and  their 

corresponding  analytic  description. 


The  next  block  of  statements  (0046  -  0108)  identifies  the  geometrical 
properties  of  the  target.  The  function  of  this  section  may  be  described 
by  consideration  of  a  particular  target.  Figure  2  shows  a  particular 
target  composed  of  a  cone-cylinder-disk.  The  geometrical  properties 
of  interest  are  the  locations  of  the  wedges,  the  wedge  angles  WA(IW), 
the  specular  directions  THSX(IW),  and  the  length  of  the  specular  line 
FLSX(IW).  This  task  is  accomplished  by  examining  the  normal  vectors 
(VNX,  VNY,  VNZ)  of  the  adjoining  surfaces  at  the  junction  between  two 
sections.  The  wedge  angle  is  then  obtained  from  the  scalar  product  of 
the  normal  vectors,  the  specular  direction  is  determined  by  checking 
for  parallel  normals  at  the  ends  of  each  section,  and  the  length  of  the 
specular  line  is  obtained  using  the  law  of  cosines  . 

The  next  block  of  statements  (0109  -  0110)  remove  the  redundancy 
in  specular  angle  which  can  occur  when  the  z-axis  is  a  specular  direction. 
For  example,  in  the  case  of  a  closed  cylinder  the  directions  0=0°,  180° 
are  specular  directions.  In  a  subsequent  test  for  a  specular  region  in 
(0190  -  0196)  the  sine  of  the  specular  angle  is  taken.  As  the  sine  is 
periodic  in  180°  care  must  be  taken  that  the  shadowed  specular  region 
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Fig.  2.  A  three  section  target-the  cone-cylinder-disk. 


is  not  included  in  evaluating  the  scattered  field.  This  is  accomplished 
by  setting  the  specular  angles  0°  and  180°  equal  to  360°.  These  specular 
directions  on  the  axis  may  then  be  evaluated  using  the  solution  for  the 
axial  caustic  given  in  Ref.  1,  and  located  at  (0242  -  0266)  in  this  program. 

The  statements  (0111  -  0115)  convert  the  computed  wedge  angles 
WA(IWD)  to  the  parameter  FN(IWD)  required  for  the  wedge  diffraction 
coefficient,  where  IWD  is  the  number  of  the  particular  wedge. 

Having  identified  the  required  geometrical  properties  of  the  target 
the  scattered  field  must  now  be  computed.  Thus  the  iteration  for  the 
incidence  angle  THT  is  set  up  (0116  -  0117)  and  the  DO  loop  entered 
(0118).  The  incidence  angle  is  incremented  and  the  propagation  vector 
of  the  incident  wave  is  computed  (0119  -  0124).  Next,  the  total  field  is 
initialized.  Then  a  DO  loop  is  entered  (0128)  in  which  the  total  scattered 
field  is  to  be  obtained  by  summation  of  the  individual  contributions  . 

First  the  location  of  the  wedge  ID  is  computed  (0129  -  0135).  Next  the 
test  parameter  BESAR  for  the  axial  solution  region  is  computed.  At  this 
point  a  transfer  to  the  solution  for  a  tip  scattered  field  is  made  if  ID  = 
NWD  (0137).  A  determination  is  now  made  whether  the  points  on  the 
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wedge  ID  at  <J)  =  0°,  180°  are  illuminated  or  shadowed  by  calling  the  sub¬ 
routine  WILLY  (0138  -  0143).  The  COMMON  block  is  now  reset  (0144)  , 
having  been  used  in  the  previous  test,  and  the  angle  between  the  incident 
ray  and  the  "side”  of  the  wedge  ID  is  computed  (0145  -  0177).  The  phase 
of  the  backscattered  field  of  the  4>  =  0°,  180°  points  on  the  wedge  ID  is 
calculated  by  the  subroutine  PHASE.  Having  completed  these  calculations 
we  now  move  to  the  computation  of  the  scattered  field  due  to  the  "wedge” 
ID. 


The  first  step  in  the  field  calculation  is  a  test  (0185  -  0189)  to  de¬ 
termine  if  the  incident  angle  THT  is  in  a  region  for  which  the  axial  solu¬ 
tion  must  be  used.  This  region  is  bounded  by  the  angle  corresponding 
to  the  first  null  of  the  first  order  Bessel  function  Jj(2ka  sinO).  If  so  a 
transfer  is  made  to  the  axial  solution.  Next  (0190  -  0197)  a  test  is  made 
to  determine  whether  the  incidence  angle  is  within  a  region  about  a  specular 
direction  specified  by  the  first  null  of  the  sin(x)/x  function.  If  so  a  trans¬ 
fer  to  the  specular  solution1  is  made.  It  is  noted  that  as  this  test  is  made 
after  the  test  for  the  axial  region  the  axial  solution  has  precedence.  At 
this  point  a  transfer  to  the  solution  for  a  tip  scattered  field  is  made  if 
ID  =  1  .  Locating  this  provision  after  the  axial  and  specular  region  test 
insures  that  the  specular  region  THSX(l)  indicated  in  Fig.  2  will  be 
included . 

Having  tested  for  regions  where  a  special  solution  is  required, 
the  scattered  field  due  to  the  wedge  discontinuity  is  now  computed 
(0201  -  0210),  and  added  to  obtain  the  total  field.  A  transfer  to  the  end 
of  the  DO  loop  is  then  made. 

The  solution  for  the  specular  field  using  the  sin(x)/x  formulation 
given  in  Ref.  2,  is  computed  in  0211  -  0238.  Again  this  field  is  added 
to  obtain  the  total  field,  and  upon  completion  transfers  to  the  end  of  the 
DO  loop. 

The  solution  for  the  axially  scattered  fields  given  in  Ref.  2  is  com¬ 
puted  in  0239  -  0272.  Again  this  contribution  is  added  to  obtain  the  total 
field,  and  a  transfer  is  made  to  the  end  of  the  DO  loop. 

The  tip  scattered  fields  are  computed  in  0273  -  0297  using  the 
physical  optics  solution  for  the  field  scattered  by  an  infinite  conducting 
cone.  Again  this  contribution  is  added  to  obtain  the  total  scattered  field. 

Statement  301  (0297)  is  the  end  of  the  range  of  the  DO  loop  over 
the  number  of  wedges.  Thus  at  this  point  the  contribution  of  the  wedge 
(ID)  to  the  total  field  has  been  computed  and  added  to  the  scattered  field. 
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At  the  termination  of  the  DO  loop  the  total  backscattered  field  EDTOT 
due  to  all  contributions  has  been  computed-  Next  the  backscattered 
field  and  the  echo  area  for  the  incident  angle  THT  are  written  (0298  -  0305), 
After  the  loop  on  incidence  angle  is  completed  the  program  is  terminated 
(0306  -  0308). 

SUBROUTINES  AND  FUNCTIONS 

Subroutine  F N ORM ( F N V X ,  FNVY,  FNVZ,  R,  THT,  PHI) 

0309  -  0331 


This  subroutine  computes  the  normal  vector  (FNVX,  FNVY,  FNVZ) 
to  the  surface  for  the  section  whose  description  is  in  the  unlabelled 
COMMON  region.  The  spherical  coordinates  of  the  point  at  which  the 
normal  is  computed  are  R,  THT  and  PHI. 

Subroutine  FCOMM(I)  0332  -  0343 


This  subroutine  shifts  the  constants  describing  the  section  I  from 
the  /DATA/ COMMON  into  unlabelled  COMMON. 

Function  RAD(THT)  0344  -  0363 


This  function  calculates  the  distance  from  the  origin  to  the 
surface  at  the  angle  THT.  The  constants  describing  the  surface  are 
stored  in  unlabelled  COMMON. 


Subroutine  WILL  Y(FWI,  THTI,  PHI,  THIN,  PHIN) 


This  subroutine  determines  whether  the  incident  wave  specified 
by  THIN,  PHIN  =  0*  tt  illuminates  the  wedge  at  location  THTI,  PHI.  If 
the  wedge  is  illuminated  FWI  is  set  equal  to  1.,  if  shadowed  FWI  is  set 
equal  to  zero  . 


Subroutine  CROSS(X,  Y,  Z,  Al,  A2,  A3,  Bl,  B2,  B3) 
0384  -  0390 


This  subroutine  computes  the  cross  product  AXB  and  returns  the 
answer  in  (X,  Y,  Z) . 
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Function  PHASE(THTI,  PHII,  THTB,  PHIB,  RB,  FK) 
0391  -  0400 


This  function  calculates  the  phase  of  the  backscattered  field  for 
an  incident  (and  scattering)  direction  given  by  THTI,  PHII  =  0  . ,  tt  ,  for 
the  wedge  location  RB,  THTB,  PHIB  =  0  . ,  tt  ,  with  a  propagation  constant 


FK. 


Function  SINXX(Y)  0401  -  0408 


This  function  computes  sin(y)/y. 


Function  DIFFl(FN,  PHI,  FK,  BETA)  0409  -  0439 


This  function  calculates  the  plane  wave  diffraction  coefficient 
given  in  Refs.  1  and  2.  FN  is  the  wedge  parameter,  PHI  the  angular 
argument,  FK  the  propagation  constant,  and  BETA  the  half  angle  of  the 
diffraction  cone. 


Function  BESLO(X)  0440  -  0465 


This  function  computes  the  zero  order  Bessel  function  JQ(x). 


Function  BESEl(X)  0466  -  0491 


This  function  computes  the  first  order  Bessel  function  Ji(x). 


Function  BESL2(X)  0492  -  0501 


This  function  computes  the  second  order  Bessel  function  Ja(x). 


INPUT  DATA 

A  typical  set  of  data  cards  is  shown  in  Appendix  II.  The  order  of 
the  cards  is  as  follows: 

Card  #1  This  card  specifies  in  format  (H5)  the  number  of  cases  to  be 
run. 
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Card  #2  This  card  specifies  the  number  of  sections  of  the  target, 
in  format  (115) . 

Card  #3  This  card  specifies,  in  order,  the  wavelength  of  the  incident 
radiation,  the  starting  incidence  angle,  the  final  incidence 
angle,  the  increment  in  angle,  and  the  polarization  factor. 

The  wavelength  is  in  meters,  the  angles  in  degrees,  and  the 
polarization  factor  is  (  +  1.)  for  E^  polarization  or  (-1.)  for  E(J) 
polarization.  The  format  is  (5F10.5). 

Next  N  cards  (N=number  of  target  sections).  These  cards  specify  in 
order,  the  initial  and  final  angular  boundaries  of  the  section 
(in  radians),  and  the  constants  Aj  .  .  .  .  .  The  format  is 

(8F10.5). 

If  the  number  on  card  number  1  is  different  from  1,  the  cards  #2  to 
(#3  +  N)  are  repeated.  This  is  necessary  if  any  of  the  following  options 
are  desired. 

1.  A  change  in  incidence  angle  range. 

2.  A  change  in  wavelength. 

3.  A  change  in  polarization. 

4.  A  change  in  target. 

It  is  also  possible  to  re-write  the  control  statements  in  the  pro¬ 
gram  to  provide  more  versatility  in  the  control,  that  is,  to  provide  for 
automatic  looping  over  wavelength,  or  polarization.  The  form  presented 
here  is  convenient  as  most  experimental  data  is  taken  in  the  form  of  an 
echo  area  pattern,  rather  than  as  echo  area  versus  frequency.  A  flow 
diagram  of  the  program  is  given  in  Appendix  IV . 


CONCLUSIONS 

This  report  is  not  intended  to  completely  specify  each  detail  of 
the  computer  program.  It  is  intended  to  provide  information  about  the 
program  which  is  necessary  to  the  use  or  modification  of  the  program. 
The  accuracy  of  the  program  has  been  tested  for  the  case  of  E^  (parallel) 
polarization  as  reported  in  Ref.  1.  Tests  for  E<J>  (perpendicular)  polari¬ 
zation  have  not  been  completed,  but  some  results  are  given  in 
Appendix  III . 
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APPENDIX  I 
COMPUTER  PROGRAM 
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$execute 
$i8 job 

$I8FTr  WEDGD 
COMPLEX 


I  8  JOB 

GOiMAP 

L I  ST  «  NODECK 

EDTOT  t PHASE • PH WU « PHWL « WDU«  WDL • EDU  «  EDL %  ESX « EBESS  <  ET IP 


3 

4 

5 

6 
7 

8 
9 
1  0 
1  1 
12 
I  3 
1  4 
1  5 
1  6 

17 

18 
19 


COMPLEX  D I FF I  i WBS «  P3ESS 
COMPLEX  WBSP«WBSM 

COMMON  RAI  »RA3iR81i R A9  «  P A  I  0  <  R  A  1  1 /OAT A/ARl  ( 20 )  « AR3 ( 20  >  «BR 1  ( 20  >  « AR9 ( 
C20) «  AR 1 0 ( 20 ) «ARI 1 (20) 

D I  MENS  ION  THSX ( 20 )  « FLSX ( 20 )  i  THT I  ( 20 )  «  THTF ( 20 ) 

D I  MENS  I  ON  VNX ( 20 )  «  VNY ( 20 )  «  VN2 ( 20 ) 

D I  MENS  ION  WA ( 21  ) »FN (21  ) 

FORMAT (  I  I  5 ) 

FORMAT ( 5F 1 0*  5 ) 

FORMAT ( 8F 1 0 • 5 ) 

FORMAT (3F15.8 ) 

FORMAT (5H  RW 1 =F I  5  •  8  t  5H  RW2=FI5.8) 

F ORM AT  ( 6H  FLSX  =  F15.8i6H  TH5X=FI5.8) 

F  ORM  AT  (  5H  WDU=2F  15*8*  5H  WDL  =  2F  1  5  •  8  «  5H  RDF=FI5.8) 

FORMAT ( 4F 1 5*8 ) 

FORMAT (  l F I  5.8 ) 

FORMAT ( 6F 1  5  •  8 ) 

FORMAT ( 6H  THT D =F 1 5 . 8 • 5H 
THTD=F15.8*7H 
ST T I =F I  5.8  ) 

ARSX  =  F 1 5 • 8 ) 

PSIU=F15.8,6H 
FWIU=F15.8*6H 
THT  W  =  F 1 5 • Q • 4H 


FORMAT (6H 
FORMAT (6H 
FORMAT ( 6H 
FORMAT ( 6H 
FORMAT ( 6H 
FORMAT ( 6H 
FORMAT ( 27H 
FORMAT ( 22H 
READ ( 5« 1 ) 


ESX  =  2F 1 5 • 8 ) 
ESESS=2F15.8 ) 


PS  I L - F 1 5 • 8 ) 

FWIL=F I 5  •  8 ) 

RW  =F15*8«4H  AW  =  F 1 5*8 i 4H  ZW=F15«8) 
PERPENDICULAR  POLARIZATION) 

PARALLEL  POLARIZATION) 

NCS 


DO  600  NCS I = 1 «NCS« 1 


READ ( 5  %  1  ) N 

READ (5*2 )  W  A  VE  «  THTS  * THTE iDTHT» POL 
WRI TE ( 6«  2  )WAVE • THTS « THTE  «  DTHT  %  POL 
I F ( POL .LT  .0,  )  WR I TE ( 6«  18 ) 

I F ( POL .GT .0 . )  WRI TE ( 6« 19 ) 

READ(5*3)  ( THT I  (  1RD)  «  THTF (  IRD)  «ARl  (  1 RD )  «  AR3(  I RD )  « 8R  I  (  I RD )  •  AR9  (  I RD )  « 
CARI 0( I RD ) • ARl I (IRD) • I RD= 1 «N ) 

WR  I  TE  (  6*  3  )  (  THT  I  (  J  )  •  THTF  (J)  «  AR  1  ( J  )  «  AR3 ( J )  « t3Rl  (J)  i  AW9  (  J  )  «  ARj  O  (  J  )  i  ARl 
Cl ( J) • J=1 «N) 

PI =3. 1415927 
TP=2.*PI 


PI2=PI/2. 

F<  =  TP/ WAVE 
DEGR A  D=  0*0 1 745329 
C  WEDGE  DETERMINATION 

C  NUM8ER  OF  TARGET  SECTIONS=N 

C  CLOSED  TARGET  ASSUMED 

NWD  =  N+ I 

dlarge=o. 

D I  A  =  0 • 

TC T  =0  *90 
DO  20 1  I W= 1  «  N  « I 
THSX (  I W  )  =  TP 
FLSX (  I W  )  =0 • 

THT 1 =THT I ( I W ) 

T HT 2=  THTF (  I W  ) 

CALL  FCOMM (  I W  ) 

RW 1 =RA  D ( THT 1  ) 

RW2=RAD ( THT2 ) 

CALL  FNORM (FXI  ,FY  1  «  FZ  1  «  RWl  «  THT  I  «0.  ) 
CALL  FNORM (FX2  •  F  Y2  <  FZ2  •  RW2 «THT2«0«  ) 
I F ( I W .FQ. 1 )  GO  TO  202 


0000 
0001 
0002 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
0010 
001  1 
0012 
001  3 
0014 
001  5 
0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
004  I 
0042 
0043 
0044 
0045 
0046 
0047 
0048 
0049 
0050 
0051 
0052 
0053 
0054 
0055 
0056 
0057 
0058 
0059 
0060 
0061 
0062 
0063 
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1W1=!W-1  0064 

T  HT  3  =  T HTF (  I W I  )  0065 

CALL  FCOMM(lWl)  0066 

RW3=RAD ( THT3 )  0067 

CALL  FNORM(FX3.FY3*FZ3.RW3.THT3.3. )  0068 

SPRD=FX1 *FX3+FYl*FY3+FZl *FZ3  0069 

V NX (  I W  )  =FX 3  0070 

V NY (  IW  )  =FY3  0071 

VNZ (  1 W  )  =FZ3  0072 

DANG=ARCOS ( ABS (SPRD ) )  0073 

I F ( SPRD • LT • 0 . )  DANG=P1-DANG  0074 

W A ( I W ) =P 1 -DANG  0075 

IF ( I W.EQ.N >GO  TO  2C3  0076 

207  CONTINUE  0077 

C  CHECK  FOR  SIN(X)/X  REGION  0078 

SPX  =  FX 1  *FX2+FY 1 *F Y2  +  FZ 1 *FZ2  0079 

SPX  =  ABS ( SPX  )  0080 

IF  (SPX.GT.TCT )GO  TO  204  0081 

GO  TO  500  0082 

202  CONTINUE  0083 

C  CHECK  FOR  '//EDGE  ANGLE  ON  THE  END  0084 

V NX (  I  )  =FX 1  0085 

VNY  (  1  )  =  F  Y 1  0086 

VNZ (  1  )  =  FZ 1  0087 

W A (  1  )=P1-2.*ATAN2(FX1  *FZ  1  )  0088 

GO  TO  207  0089 

203  CONTINUE  0090 

C  CHECK  FOR  WEDGE  ANGLE  ON  THE  END  0091 

FZ2A=ABS(FZ2 )  0092 

VNX ( NWD ) =FX2  0093 

VNY(NWD)=FY2  0094 

VNZ (NWD ) =FZ2  0095 

WA(NWD)=P1-2.*ATAN2  (FX2.FZ2A  )  0096 

GO  TO  207  0097 

204  CONTINUE  0098 

THSX (  IW  )  =  AT AN2 (FX1  iFZ 1  )  0099 

THHX  =  THSX (  1 W )  0100 

I  F (  I  W • GT • 1  .AND*  ABS(THHX ) .EQ.O.  )  THSX(1W)=P1  0101 

FLSX  (  1  W  )  =  SORT  (  ABS  (  RvV  1  *R'W1  +W  W2*« W2~2  •  *» W  1  *RW2*C0S  (  THTp-THT  1  )  )  )  0  102 

GLSX=FLSX ( 1 W )  0103 

WRITE (6.5)  RW1.RW2  0104 

WR I TE (6*5)  THTI.THT2  0105 

WR I TE ( 6 . 6 )  GLSX.THHX  0106 

500  CONTINUE  0107 

201  CONTINUE  0108 

IF ( THSX  (  1  )  .LT  .0. 1  )  THSX(1)=TP  0109 

1F(A6S(THSX(N)-P1).LT.0.1)THSX(N)=TP  0110 

DO  208  I WD“ 1  *  NWD  *  1  0111 

FN ( I WD ) = ( TP-WA ( 1WD ) ) /P 1  0112 

FNWR=FN ( 1WD )  0113 

W  R 1 T  E (6*9 ) F  N  W  R  01  14 

208  CONT 1NUE  0115 

FLOOP  =  ABS (  ( THTS-THTE ) /DTHT  )  0  116 

LOOP =FL OOP +1 •  0117 

DO  400  LP= 1  *  LOOP .  1  0118 

FLP=LP  0119 

THT  D= (FLP-1  •  )  *DTHT  0  120 

THT=DEGRAD*THTD  0121 

VX=-SIN(THT)  0122 

VY  =  0 »  0123 

VZ=— COS (THT )  0124 

C  PARALLEL  POLARIZATION  0125 

EDTOT  = (O. . J.  )  0126 

1 DSX  =  0  0127 
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DO  301  1 D  = 1  *N'WD*  1  0128 

THTW=THT I ( ID)  0129 

CALL  FCOMM (ID)  0130 

IF ( 1D.EQ.NWD)  THTW=THTF(N)  0131 

RW=RAD ( THTW )  0132 

AW=RW*S1N(THTW)  0133 

ZW=RW*COS ( THTW )  0134 

WR 1 TE ( 6 1 17)  THTW*RW* AW*ZW  0135 

0ESAR  =  2  «*FK*AW  0136 

IF ( 1D.EQ.NWD )  GO  TO  305  0137 

I F ( THT .LT  .P 1 2 • AND.  1 D*GT •  1  )  CALL  FCOMM(lD-l)  0  138 

CALL  'W  I  LL  Y  (  F  W  1  ♦ THTW  iO#  ♦ THT *0*  )  0  139 

FW1U  =  F’WI  0140 

CALL  WILLY  (  F  W  1  1  THT'a'  iPl  »  THT  <0*  )  014  1 

FW1L=FW1  0142 

WR 1  T  E  (  6 «  16)  FWIU^FWIL  0143 

CALL  FCOMM( ID)  0144 

306  CONTINUE  0145 

FIND  INCIDENCE  ANGLE  ON  WEDGE  0146 

INCIDENCE  VECTOR  VX*VY*VZ  0147 

VXM=-VX  0148 

V  Y  M  =  —  V  Y  0149 

V  Z  M  =  -  V  Z  0180 

1F( THT.GT.PI2  )  GO  TO  309  0151 

SPRD= VXY* VNX ( I D ) +VZM*VNZ ( ID)  0152 

CALL  CROSS  (BX*BY«BZ»VXM,\/YM*VZM«VNX(  1  D  )  ♦  VNY  (  ID)  « VNZ  (  ID)  )  0  153 

DANG= ARC OS ( SPRD )  0154 

S 1 GN  =  - 1  .  0155 

1 F ( BY • LT • 0 • )  S1GN=1.  0156 

PS  I U  =  P 12  +  S 1GN*DANG  0 1 57 

SPRD=-VXM*VNX (ID) +VZM*VNZ (ID)  0158 

CALL  CROSS  U3X  *8 Y  *8Z ♦ VXM, VYNU  VZM,-VNX (  ID  >  *VNY  <  ID  >♦ VNZ  <  ID )  >  0159 

DANG=ARCOS (SPRD)  0160 

S 1 GN - 1  *  0161 

1F(BY.LT.C*)  SlGNr-i.  0162 

PS  1 L  =  P 12+S 1 GN*D ANG  0 1 63 

GO  TO  310  0164 

309  CONTINUE  0165 

CALL  FNORM(FTX,FTY«FTZ*RW* THT  '  «  0*  )  0166 

SPRD=VXM*FTX+VYM*FTY+VZM*FTZ  0 1 67 

CALL  CROSS  (  BX  »RY  1  HZ  1  VX’^li  VYM  »  VZM  »FTX  1  FT  Y  iFTZ  )  0168 

DANG=ARCOS ( SPRD )  0169 

S 1 GN= 1 •  0170 

1 F ( B Y • L  T • 0 •  )  S 1 GN  =  -  1  •  0171 

PS  I U  =  P I 2+S IGN*DANG  0  I  72 

SPRD=-VXM*FTX+VYM*FTY+VZM*FTZ  0 1 73 

DANG  =  ARC OS ( SPRD  )  0  174 

PS  I L=P I  2 -DANG  0  175 

310  CONTINUE  0176 

WRITE (6« 15)PS1U»PSIL  0177 

C  DETERMINE  PHASE  of  WEDGE  0178 

PHWU  =  PHASE  (THT  «0i  «  T  H  T  «  0  •  »RW  «FK  )  0  179 

PHWL  =  PHASE ( THT ♦  0  •  •  THT  A  «  P 1  «  RW  *  FK  )  0180 

P  SUP  =  2 • *PS 1 U  0181 

PSUM=0.  0 1 52 

PSLP=2.*PSIL  0183 

PSLM=0.  0184 

C  8ESSL0  TEST  0185 

ST TO =3. 831 7/BESAR  0186 

STT 1 =S IN ( THT )  0137 

WR I TE  ( 6  *  13)  STT 1  0188 

IF  (STT I .LT .SfTC ) GO  TO  302  0189 

c  S I N ( X ) / X  TEST  0190 

1 D  X  = ID  0  191 
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T 1 =  ASS ( THSX (  1 DX  >-THT )  0  192 

I F ( T 1  *  GT • P 1 2 )  T 1 =  P 1 2  0193 

ARSX  =2  .*FK*FLSX(  1DX  )*SIN(T 1  )  0194 

WR1 TE ( 6*  1 4  )  ARSX  0195 

1F( ARSX.LT.Pl )G0T0  304  0196 

303  CONTINUE  0197 

NCHECK  =  303  0198 

WR 1 TE (6*1)  NCHECK  0199 

lF(lD.EO.l)  GO  TO  305  0200 

C  CALCULATE  DIFFRACTED  FIELDS  0201 

WDU  =  D1FF1  ( FN (  ID)  »PSUPiFK,Pl2  )+POL*DlFFl  (FN( ID)  «PSUM,FK*P12  )  0202 

IF ( 1D.EO. 1DSX )  WDU= ( 0 • • C • )  0203 

WDL=DIFFl  (FN(  ID)  *  PSLP  *  FK  *  P  I  2  )+POL*DlFFl  (FN(  ID)  <  PSLM  »FK»  P  12  )  0204 

RDF  =  SQRT  (At3S(  AW/(2**SIN(THT  )  )  )  )  0205 

WR  1  TE  (6*7)  WDU  ♦  WDL  ♦  RDF  0206 

EDU=RDF*PHWU*WDU*FW  lu  0207 

EDL=RDF*PHWL*WDL*FW I L  0208 

EDTOT=EDTOT+EDU+EDL  0209 

GO  TO  3^1  0210 

304  CONT I NUE  021 1 

NCHECK =304  0212 

WR 1 TE ( 6  *  1  )  NCHFCK  0213 

COMPUTE  S I N ( X  )  / X  CONTRIBUTION  0214 

COMPUTE  PEAK  RETURN  USING  AVERAGE  OF  GAUSSIAN  CURVATURE  0215 

R W 1 =  RW  02 1 6 

RW2=RAD ( THTF ( ID) )  0217 

A  W 1 =  A  W  0218 

AW2=RW2*SIN (THTF( I D ) )  0219 

TAW=A5S ( AW1 -AW2 )  0220 

1F(TAW*LT*0«1  )  GO  TO  307  0221 

RHI =AW 1 /S IN ( THT )  0222 

RH2=AW2/S IN( THT )  0223 

RH 1 =ABS ( RH 1 )  0224 

RH2=ABS(RH2)  0225 

RSH1 =SORT (RHI )  0226 

RSH2  =  S0RT (RH2  )  0227 

PEAK= ABS (RSH 1 *RH 1 -RSH2*RH2 ) *FLSX ( 1 D ) *1 . 4 1 4/ ( 3«*ABS (RHI -RH2 ) )  0228 

PEAK=FK*PEAK  0229 

GO  TO  308  0230 

30  7  PEAK=S0WT ( AW  1 /2*  ) *FK*FLSX (ID)  0231 

308  CONTINUE  0232 

ESX=-PEAK*S I NXX ( ARSX ) /SORT ( TP*FK )  0233 

ESX  =  ESX*PHASC ( THT • 0 • • THT W «  0 •  «  RW ♦ FK )  0234 

1 DSX  = 1 D  + I  0235 

WR I TE (6*11) THTD.ESX  0236 

EDTOT=EDTOT+FSX  0237 

GO  TO  301  0238 

302  CONTINUE  0239 

NCHECK  =302  0240 

WRITE (6*1)  NCHECK  0241 

C  COMPUTE  FIELD  USING  BESSEL  FUNCTIONS  0242 

S8SS=0.5  0243 

PKSPC=FK*A W* AW  0244 

PKBSS=Aa  0245 

ARGBS=6ESAR*S I N ( THT )  0246 

I  F  (  Al3S  (  ARGBS  )  *GT  •  0  •  05  )  SBSS=BESL  1  (  ARGBS  ) /ARG13S  0247 

CALL  FNORM ( F 1 X , F 1 Y  *  F I  7  % RW • THT  W «  0 •  )  0248 

P  SUB  =  2 • * ( P 1 2-ARCOS ( AHS (F I Z  )  )  )  024  9 

1 F ( THT •  L  T • P I  2  )  PSUd=2 •* ( P I  2- ARCOS ( VNZ (  ID)  )  )  0250 

WBS  =  DIFF1  ( FN (  ID)  .PSU8«FK«P12  )  0251 

P8ESS  =  PHASE ( THT *  0.  . 0 . . D •  • ZW . FK )  0252 

SHBDT  =  C ABS ( W3S ) -0. 5  0253 

IF(  A3S  (SHL3DT  )  .LT.0,05)  GO  TO  32  1  0254 

WBS=WBS*SQRT ( TP*FK )  0255 
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COST=COS (THT )  0256 

WBSP  =  D IFF  I  (FN (  ID )  »  0 • *FKf P 12 ) +DIFF I  ( F  N (  ID) »  PSUB  «  FK •  P  I  2  )  0257 

WBSM  =  D IFF  I  (FN (  ID )  • 0 • » FK »P 12 ) -D IFFI  ( FN (  ID)  *  PSUB • FK  .  P I  2 )  0258 

WBSP=WBSP*SQRT (FK*TP )  0259 

WBSM=WBSM*SQRT ( TP*FK )  0260 

WBS=-COST*COST*WBSM*SBSS+WBSP* ( SBSS-BESL2 ( ARGBS ) )  0261 

I F ( POL  •  L T  •  0 •  )  WBS  =  C0ST*C0ST*WUSP*SBSS-WBSM*<SBSS-I3ESL2< ARGBS )  )  0262 

EBESS= ( 0, ,-l • ) *PKBSS*PBESS* WBS  0263 

IF<FWIU.EQ*0. .AND.FWIL#EQ*0. )  EBESS= (0«iO. )  0264 

I F ( FW I U • GT • 3 • . AND.FWIL.EQ.O.  )  EBESS  =  0 • 5*EBFSS  0265 

GO  TO  322  0266 

321  CONTINUE  0267 

EBE  SS  = ( G.  ,-i  .  )*PKSPC*PBESS*SBSS  0268 

322  CONTINUE  0269 

WR I TE ( 6  *  12)  THTD«EBE3S  0270 

EDTOT=EDTOT+EBESS  0271 

GO  TO  30 I  0272 

305  CONTINUE  0273 

NCHECK  =  30  5  0274 

W  R  I  T  E  (  6  *  I  )  NCHECK  0275 

COMPUTE  TIP  CONTRIBUTION  AT  ENDS  0276 

TO  FIRST  ORDER  TIP  CONTRIBUTIONS  ARE  NEGLEGIBLE  0277 

THTCN=THT  0278 

IF ( THT .GT .P 12  )  THTCN  =  PI-THT  0279 

HFC A  =  0 • 5*W A (  ID)  0280 

T CN A  =  P I 2-HFCA  028  I 

T CN A A  =  ABS ( THTCN-TCNA  )  0282 

IF ( TCNAA.LT .0. I  745)  GO  TO  312  0283 

I F ( ABS ( HFCA-P 12 ) .LT  .0.  I  745  )  GO  TO  312  0284 

T  AN A  =  T  AN (HFC A  )  0285 

T  ANA2  =  T  AN  A*T  ANA  0286 

TANTH=TAN ( THTCN )  0287 

TANT2  =  TANTH*T  ANTH  0288 

COSTH  =  COS ( THCN  )  0289 

C0STH3=C0STH*C0STH*C0STH  0290 

ETIP  =  -WAVE*TANA2/(8.*PI*C0STH3* ( I  . -T ANA2*T ANT2  )  )  0291 

E  T'l  P  =  ET  I  P*PH  ASE  <  THT  *0.  t  0.  t  ZtffFK  )  0292 

IF (POL.LT. J#  )  E  T I P  =  — ET I P  0293 

GO  TO  31  I  0294 

3  12  ET IP= ( 0.  %  0.  )  0295 

31  I  CONT INUE  0296 

EDTOT=EDTOT+ETI P  0297 

301  CONTINUE  0298 

C  OUTPUT  SECTION  0299 

WR I TE ( 6  *  4  )  THTD*EDTOT  0300 

E  M AG  =  C ABS (EDTOT )  030  1 

S IGMA=2«*TP*EMAG*EMAG/ ( WAVE* WAVE )  0302 

S I GM AL  = I  0 . * ALOG I  0 ( SI GMA  )  0303 

WRITE (6% 4  )  S I GMA . S I GM AL  0304 

400  CONTINUE  0305 

600  CONTINUE  0306 

STOP  0307 

END  0308 

SIBFTr  FFNRM  L I  ST  *  NODECK  0309 

SUBROUTINE  FNORM(FNVX.FNVY,FNVZ.R.THT.PHI )  0310 

C  INPUT  R»THTiPHI  0311 

COMMON  API  . AR3. B« 1  *  AR9. API  0 . AR I  I  0312 

ST  =  S I N ( T  HT )  0313 

CT  =  COS ( THT  )  0314 

S P  =  S I N ( PH  I  )  0315 

CP=COS (PHI)  0316 

U  =  AR I *ST*ST  +  AR3*CT*CT  +  6RI *ST*CT  0317 

V  =  AR9*CT  +  AR I  0  *S T  0318 

UTH=2 • * ( AR I -AR3 )*ST*CT+BRI * ( CT*CT-ST*ST )  0319 
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VTH=~AR9*5T+ AR l 0*CT  0320 

F1=2.*R*U+V  0321 

F2=R*UTH+VTH  0322 

F  X  =  ST  *CP*F 1 +CT*CP*F2  0323 

FY=ST*SP*F1 +CT*SP*F2  0324 

FZ=CT*F 1-ST*F2  032b 

F N  =  SORT ( F X*F X  +F  Y*F Y  +F  Z*F Z )  0326 

FNVX=FX/FN  0327 

FNVV=FV/FN  0328 

FNVZ=FZ/FN  0329 

RETURN  0330 

END  0331 

MBFTC  F  COMM  •  LIST. NODECK  0332 

SUBROUTINE  FCOMM(I)  0333 

COMMON  RA 1  .  RA3  ,  RBI  .  RA9 % RA I  0 . RA 1  1 /DA T A/AR 1  (20  )  * AR3 ( 20 )  , 3R 1  ( 20 ) ♦ AR9 (  0334 

C2  0 )  . AR 1 C ( 2 J  )  .  AR1  1  ( 2  0  )  0335 

RA 1 = AR 1 (I )  0336 

R A  3  =  AR3 (  I  )  0337 

RBI  =  BR 1  (  I  )  0338 

R  A  9  =  A  R  9 (  I  )  0339 

R A  I  0  =  AR 1 0  (  I  )  0340 

R  A  1  I =  AR I  1(1)  0341 

RETURN  0342 

END  0343 

S  I BF  T C  RADD  LIST. NODECK  0344 

FUNCTION  RAD(THT)  0345 

C  CALCULATE  R  0346 

COMMON  AR 1  * AR3.BRI  . AR9*ARl OVARII  0347 

R  =  0  *  0348 

ST  =  S I N ( THT  )  0349 

C  T  =  COS ( THT  )  03b0 

C 1  =  API *ST*ST  +  AR3*CT*CT  +  bP1 *ST*CT  03b  I 

C2=AR9*CT+ARlO*ST  03b2 

C3= AR1 1  03b3 

IF (Cl • EO • 0 • ) GO  TO  II  03b4 

ARG=SQRT ( C2*C2-4 .*C 1 *C3 )  0355 

R= ( -C2+ARG ) / ( 2.*C 1 )  0356 

R2= ( -C2-ARG )/ ( 2 • *C 1 )  0357 

I F ( R2.LT .R. AND.R2.GT .0. )  R=R2  0358 

GO  TO  12  0359 

1  1  R  =- C3/C2  0360 

12  RAD=ABS(R)  0361 

RETURN  0362 

END  0363 

MSFTr  WILLY.  LIST. NODECK  0364 

SUBROUT  I NE  W I L L Y ( F W I  . T H T I  .PHI  .THIN.PHIN)  0365 

C  INPUT  THT 1  . PH  I  . TH I N . PH  I N  0366 

C  F W I =0 •  WEDGE  IN  SHADOW  0367 

C  F  W  I  =  1  •  WEDGE  ILLUMINATED  0368 

COMMON  API  . AR3 ♦ 0R I  . AR9 . A  R I  C  .  AR  1  1  0  369 

P I =  3 •  1 4 1 5927  0370 

P I  2= I  .5707963  0371 

R=R AD ( THT I )  0372 

CALL  FNORM (FNX.FNY .FNZ.R. THT I .PHI )  0373 

TH=THIN  0374 

X  I N  =  S I N ( TH ) *COS (PH  I N )  0375 

Y I N  =  C  *  0376 

Z  I N  =  CO  S ( TH )  0377 

SPRD  =  F  NX*X I N  +  FNY*Y I N+FNZ*Z IN  0378 

FWl = 1 .  0379 

I F ( SPRD.LT • 0 . )  FW I =0.  0380 

CONTINUE  0361 

RETURN  0382 

END  0383 


17 


o  o 


SIBFTC  CROSS.  LIST. NODECK  0384 

SUBROUTINE  CROSS (  X  .  Y  •  Z  .  A  I .A2.A3.B1. B2 .03)  0385 

X=A2*B3-A3*B2  0386 

Y=A3*B1-AI *B3  0387 

Z=A 1*B2-A2*0 I  0388 

RETURN  0389 

END  0390 

SIBFTC  PHAS.  LIST. NODECK  0391 

COMPLEX  FUNCTION  PH A SE < THT I  . PH  I  I  . TH TB . PH  I B . RB . FK >  0392 

FS= I •  0393 

TE  S  T  =  AB  S  <  PH  I  I —PH  I  6 )  0394 

IF(TEST.EQ.0. )  FS=- I •  0395 

FL  =  2.*R3*COS ( THT I +FS*THTB  )  0396 

FKL=FK*FL  0397 

PHASE  =  CMPLX ( COS (FKL  )  .S I N (FKL  )  )  0398 

RETURN  0399 

END  0400 

SIBFTr  SINXX.  LIST. NODECK  0401 

FUNCT I ONSINXX ( Y )  0402 

C  S I N ( X ) /X  FUNCTION  0403 

Z=Y*Y  0404 

SI NXX= (  I  •- (Z/6.  >* (  I  •-< Z/20 •  ) *  <  I  .- < Z/42 •  ) * < I • - ( Z/72 •  > * ( I  • - ( Z/ I  I  0 •  > *  0405 

2(  I  (Z/I 56.  )* (  I  •  - ( Z/I 82.  )  )  )  )  )  )  )  )  0406 

RETURN  0407 

END  0408 

SIBFTC  DIFCN  LIST. NODECK  0409 

COMPLEX  FUNCTION  D I FF 1  ( FN . PH  I  . FK . BET A )  04  10 

SINGLE  DIFFRACTION  FUNCTION  D I FF I  0411 

FN-WEDGE  ANGLE.  PH  I c ARGUMENT .  BETA  =  CONE  ANCLE  04  12 

COMPLEX  EXPH  0413 

PI  =3.  I  415927  0414 

TP=2.*PI  0415 

P I  4  =  P I /4 •  0416 

PIF=PI/FN  0417 

PHN=PH I /FN  0418 

PHI A=ABS(PHI )  0419 

EXPH  =  CEXP (  ( 0 •  .-1 •  ) * P I  4 )  042  0 

S I NB  =  S I N ( BETA  )  042  I 

S I NB  =  ABS ( S I N9  )  0422 

IF (S1N3.LT. .001 )  GO  TO  3  0423 

I F ( ABS (PH  I A-P I  )  ,LT. 0. I  )  GO  TO  I  0424 

T 1 =  COS ( P I F  )  0425 

T2  =  COS ( PHN )  0426 

T3=SIN(PIF)  0427 

T4=TI-T2  0428 

T  4  A  =  ABS ( T4 )  0429 

IF ( T4A.LT. .00  I  )  GO  TO  3  0430 

DIFFI =EXPH*T3/(FN*SIN3*T4 )  0431 

DIFFI  =D I FF I /SORT ( TP*FK  )  0432 

GO  TO  2  0433 

1  DIFFI =0.5/SlNB  0434 

GO  TO  2  0435 

3  D I FF I = I COOO.  0436 

2  CONTINUE  0437 

RETURN  0438 

END  0439 

SIBFTC  BESLO.  LIST. NODECK  0440 

FUNCTION  BESLO(X)  0441 

C  (X.GT.O. )  0442 

Z  =  ABS  ( X  )  0443 

Y  =  X / 3 •  0444 

Y2  =  Y  * Y  0445 

Y  3  =  Y  2  *Y  0446 

Y4=Y2*Y2  0447 
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Y5=Y3*Y2  0448 

Y6=Y4*Y2  0449 

Y8=Y6*Y2  0450 

Y 1 0  =  Y8  * Y  2  04b 1 

Y 1 2  =  Y I  0 * Y 2  04b2 

IF (Z.GT .3. )GO  TO  lC  0453 

8ESL0=l ,-2.2499997*Y2+ 1 • 2656208* Y4- . 3 1 63866* Y6 +  • 0444479* Y8- . 003944  0454 

24*Y1  0-.00'J24R46*Yl2  0455 

GO  TO  1 1  04b6 

10  CONTINUE  0457 

F= . 79788456- ( ,00000077/7 )- ( • U0552740/Y3 )- ( ,0000951 2 /Y3 >+(.00137237  04 b8 

2/Y4  )  -  (  •  w'O 0728 05/Y5  >  +  (  ,00  0  1  4  476/Y6  )  04  59 

T  =  X-. 785393 1 6- (  ,04 1 6639  7/Y )- ( . 0000 3954 /Y2 >+ (  .0 0  262573/ Y 3 >- ( . 00  054  I  04  60 

225/Y4)-(.OJ029333/Y5)+(,0C0l 3558/Y6 )  0461 

8ESL0=F*CO5( T ) /SORT (X )  0462 

I  I  CONT I NUE  0463 

RETURN  0464 

END  0465 

SlRFTr  BESL1.  LIST, NODECK  0466 

FUNCTION  BESLI(X)  0467 

C  POLYNOMIAL  APPROXIMATION  X.GT,0.  0468 

Y  =  X/3 •  0469 

Y  2  =  Y  *  Y  0470 

Y3=Y2*Y  0471 

Y4=Y3*Y  0472 

Y5=Y4*Y  0473 

Y6=Y5*Y  0474 

Y8=Y6*Y2  0475 

Y 1 O  =  Y  8*  Y  2  0476 

Y12=Y1C*Y2  0477 

I F ( X • GT • 3 •  >  GO  TO  lO  0478 

BESL  =  0. 5- • 86249985* Y2  +  . 21 093573 *Y4 - . 03954  289* Y6+ • 0  044  3  3  I  9*Y8-. C003  04  79 

21 761  * Y 1 0+ ,00001  1 09*  Y 1 2  048  0 

BESLl =X*8EoL  0481 

GO  TO  11  0482 

10  CONTINUE  0483 

F  =  *  7  97884  56+ (,0000  0  1 56/ Y  )+ ( • 0 1659667/ Y 2 )  +  (  .00  0 1 71 05/Y3 ) -(.00249511  0484 

2/Y4  )+ ( • 001  1 3653/Y5  > - ( .00020033/Y6 >  04  85 

T  =  X-2 • 356 1 94  4  9+ ( .  1 24  996 1 2/Y  )  + ( , GG00565  0/Y2 ) - (  .00637879/Y3 ) + ( .00074  0486 

2348/Y4  )  +  ( .000 79824 /Y5  >- ( . 000291 66/ Y6  >  0487 

8ESL 1 =F*COS ( T ) /SORT (X )  0488 

II  CONTINUE  0489 

RETURN  0490 

END  0491 

SIBFTr  BESL2.  LIST, NODECK  0492 

FUNCTION  B  E  3  L  2 (X )  0493 

C  ( X.GT . 0  •  )  0494 

IF (X.EQ.O. ) GO  TO  lC  0495 

BE  SL2  = ( 2 • *3ESL 1  ( X ) /X ) -BESLO ( X )  04  96 

GO  TO  11  0497 

10  HESL2  =  0 .  0498 

1 1  CONT 1 NUE  0499 

RETURN  0500 

END  0501 

3D A  T  A  0b02 
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APPENDIX  II 
INPUT  DATA 


00002 

00003 

1.0 

0. 

180. 

1.0 

0. 

.785398 

0. 

0. 

.785398 

2.617 

0. 

0. 

2.617 

3 . 1415927 

0. 

0  . 

00003 

1.0 

0. 

180. 

1 .0 

0. 

.785398 

0  . 

0. 

.785398 

2.617 

0. 

0. 

2.617 

3 . 1415927 

0. 

0. 

-1 . 

0.  2.0 

0.  0. 

0.  -1.0 

1.0 

0.  2.0 

0.  0. 

0.  -1.0 


1.0  -2.0 

1.0  -1.0 

0  -2.0 


1.0  -2.0 
1.0  -1.0 

0.  -2.0 


20 


APPENDIX  III 

COMPUTER  TESTS  FOR  E<j>  POLARIZATION 


A  limited  amount  of  data  has  been  obtained  for  the  case  of  E<j) 
polarization.  The  computed  results  are  compared  with  measured  data 
for  the  cone  and  double  cone  in  Figs.  3  and  4.  The  data  for  the  cone 
is  that  of  Keys  and  Primich5  and  the  double  cone  data  is  from  Eberle 
and  St-  Clair.6  For  the  case  of  the  cone  the  greatest  error  is  in 
the  near  forward  region  where  the  single  diffraction  solution  is  used. 

The  results  in  this  region  agree  with  the  single  diffraction  results 
of  Bechtel.7  The  results  in  this  region  could  be  improved  by 
extending  the  edge  current  formulation  for  the  axially  scattered  fields 
to  this  region.  In  the  case  of  the  double  cone  the  calculated  echo  area 
is  in  error  in  the  region  about  0  =  90°  .  This  is  a  region  where  the 
effect  of  edge  type  creeping  waves1  at  the  junction  should  be  considered. 
These  test  results  indicate  that  further  extension  of  the  theory  and  com¬ 
puter  program  is  necessary  to  remove  these  sources  of  error.  However, 
for  targets  which  are  large  in  terms  of  wavelength  the  accuracy  will 
improve,  as  both  the  edge  current  and  creeping  wave  effects  would 
diminis  h . 
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180 


Fig.  3.  Echo  area  of  a  cone-E^  polarization. 
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Fig*  4.  Echo  area  of  a  double  cone-Eep  polarization. 
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APPENDIX  IV 

WEDGE  DIFFRACTION  COMPUTER  PROGRAM 
FLOW  DIAGRAM 


The  flow  diagram  for  the  wedge  diffraction  Computer  program 
WEDGD  is  given  in  Fig.  5.  This  flow  diagram  has  also  been  presented 
in  Ref  •  1 . 
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Wedge  D  »  rF/?Acr/o^  Co^p uter  program 
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